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Abstract 

We review our recent proposal for a universal description of generic single-component 
viscoelastic systems with a single relaxation time. Foliation preserving diffeomor- 
phisms are introduced as an underlying symmetry which naturally interpolates be- 
tween the two extreme characters of elasticity and fluidity. The symmetry is found 
to be powerful enough to determine the dynamics in the first order of strains. 
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1 Introduction 



Viscoelasticity (see, e.g., [HE]) is a unifying notion which contains sohd and fluid as hmiting 
cases. It was defined by James Clerk Maxwell such as to describe the materials that behave 
as elastic solids at short times scales and as viscous fluids at long time scales [3], H]. 

In order to grasp the idea of this definition, we consider a material (like a chewing gum) 
consisting of many particles bonding each other. They stay at their equilibrium positions in 
the absence of stresses (as in the leftmost of Fig. [1]), and we now suppose that an external 
force is applied to deform the material. An internal stress is then produced in the body, 
and according to the definition, this internal stress can be treated as an elastic force at least 
during short intervals of time. However, if we keep the deformation much longer than some 
period of time (characteristic to each material and to be called the relaxation time), then 
the bonding structure changes to reduce the free energy, and the internal stress vanishes 
eventually. The point here is that the two configurations (the center and the rightmost of 
Fig. [T]) exhibit the same shape (i.e. same configurations for the positions of particles) but 
have different bonding structures. 

shear 
deformation 

>■ 

Figure 1: Shear deformation of a viscoelastic material 

We thus see that the bonding structure itself is dynamical for viscoelastic materials, and 
is lead to introduce a new dynamical variables (to be called intrinsic metric [5]) to specify 
the bonding structure. It was claimed in [5] that viscoelasticity may have a universal 
description if we define the intrinsic metric properly and assume the invariance of the 
system under foliation preserving diffeomorphisms (FPDs). In fact, it was shown there 
that FPDs naturally interpolate between the two extreme limits of elasticity and fluidity, 
and are powerful enough to determine the dynamics in the first order of strains. 

The main purpose of the present article is to give a brief review of [5] in a concise way, 
and also to present some additional stuffs which were not mentioned in [5]. 

Our discussions will be constrained to single-component materials. Note that, even for 
this case, there can be two relaxation times, one being for shear deformations and another 
for bulk compressions (or expansions). In this article, we set the relaxation time for bulk 
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compressions to be infinite, so that the bulk stress does not undergo relaxations (see Fig.|2]). 
This ensures that the system exhibits fluidity at long intervals of time. 




compression 
► 




no change 
► 




Figure 2: Compression of a viscoelastic material 



2 Geometrical setup 

We describe a viscoelastic material within the framework of world-volume theory with the 
world- volume coordinates ^ = (^'^) (a = 1,2,3) (labeling the material particles of a body) 
and with the target space coordinates X = (X*) (i = 1,2,3). The shape of the material 
at time t is then specified by giving the functions X*(,^, t). In order to further describe the 
bonding structure, we introduce an intrinsic metric habi^,t) as follows. 



We first recall that our viscoelastic materials are defined to be elastic at least for a very 
short time. This means that for any (sufficiently small) portion in a material, its natural 
shape can be defined at each moment as the shape that would be taken when all stresses 
were virtually removed (see Fig. H]). 

This observation enables us to define two kinds of metrics on the world-volume at each 
time t (see Fig. |5]). The first one is the induced metric habi^,t) which measures the length 
in the real shape: 




Figure 3: World- volume of a viscoelastic material 



(1) 
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Figure 4: Definition of tlie natural sliape for a small portion 



where 

habiC, t) = e^(e, t)el{C, t) (ei(e, t) = daX\^, t) : dreibein) . (2) 
Another is the intrinsic metric hab{^,t) which measures the length in the natural shape: 

ds^ = hab{^,t)dCd^\ (3) 
The difference between them defines the strain tensor, 

e,,{^,t) = ^{hab{^,t) -hab{^,t)). (4) 

In fact, for an elastic material at rest without stress, we can take the coordinates ^ = {^^) 
such that hab = Sab- Then, when the material is slightly deformed under stresses, X°-{C,, t) = 
^'^ + u'^{^, t) (u": displacements), we have hab = daX^dbX^ = 5ab + daUb + dbUa + 0{u^), so 

that Eab = {l/2){daUb + dbUa) + 0{u^). 




Figure 5: The induced and intrinsic metrics 



We denote by po the mass density in the absence of strains. In this article, we assume 
that po is a constant independent of ^ and t. Then the mass contained in a volume element 
d^^ = d^^d^'^d^^ is given by pQ^/fid'^^. The mass density p in the real three-dimensional 
space is given by the mass conservation pVhd^^ = po^/hd^^ as p = (v^/v^) po- 

Let us consider various deformations of a 2D material. Here are comoving 

coordinates attached to material particles. 
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(a) shear deformations: Suppose that a material in equilibrium is deformed with an 
external force as in Fig. |6l The two metrics hat and hab agree before the deformation. When 
the material is deformed, the induced metric hab varies, expressing the change of the shape. 
In contrast, the intrinsic metric hab retains its original form just after the deformation. 
This is because the material should exhibit elasticity for a short period of time, and the 
natural shape at that moment is its original form in equilibrium. At sufficiently later times, 
hab gradually changes its values and comes to equal the induced metric. 




Figure 6: Shear deformation with the angle 6 

(b) bulk compressions: Suppose that a material in equilibrium is compressed with a 
factor A (< 1) as in Fig. [71 The two metrics hab and hab again agree before the deformation, 
and the induced metric hab changes after the deformation. However, according to our 
assumption that the bulk compression (or expansion) does not undergo relaxation, the 
intrinsic metric hab retains its original form. 

3 Foliation preserving diffeomorphisms 

So far the world-volume coordinates have been attached to a material, labeling the 
material particles. However, we can introduce arbitrary coordinates C, at each time t, 
independently of the labeling of material particles. In this case we need to specify the 
relative motion of the material to the coordinates, but this can be carried out by introducing 
a vector field A/" = N°'{^,t)da (see Fig. E]) [3]. More precisely, the vector field is defined 
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Figure 7: Compression by the factor A (A < 1) 



such that the material particle located at ^ at time t is supposed to move to the position 
_j_ iV"(^, t)dt after the time interval dt. 




Figure 8: Definition of A^" 



With N°'{^,t), the system comes to have the gauge symmetry of foliation preserving 
diffeomorphisms (FPDs), reflecting the redundancy in the description. The set of FPDs 
consist of the (3 + l)-dimensional reparametrizations that do not move the time slices, 
and is generated by the infinitesimal transformations of the form (see [3] for references of 
mathematical details): 

SC{U) = e^{^,t), 6t{^,t) = 0. (1) 
The fields X*, hab and hab are intrinsically three-dimensional and thus transform covariantly 
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as 

5X' = e''daX\ (2) 

5hab = dae'' hcb + dbe" hac + dchab , (3) 

Shab = dae" hcb + dbe" Kc + e'' djiab , (4) 

while the transformation of the shift vector has an additional inhomogeneous term as 

57V« = -d^e"" + dcN" - t . (5) 

The time derivatives of a FPD-covariant tensor Tj^la^..'. are no longer covariant, but can be 
made so by introducing the covariant derivative [5]: 

D d 

rpb\b2- - — rpblb2- - i r rpblb2--- (c\ 

Jjj. CLia2 - — Q^-^aia2- - ~^ N a]_a2- ■ ■ ' I 

where Cj^ is the Lie derivative with respect to M = N°'da- 

For example, the velocity v^{C,,t) and the acceleration a^{C,,t) of the material particle 
located at t) are given, respectively, by 

v^ = ^=r + N-daX' , = ^ = = ^' + N^d^v' . (7) 

This can be easily confirmed by noting that the material particle at (^, t) moves after the 
time interval dt to the point + Ndt, t + dt): 

X\i + Ndt, t + dt)= X\i, t) + [X\i, t) + iV-(e, t) daX\i, t)] dt = t) + v'{^, t) dt . 

(8) 

As for the metrics hab and hab, their covariant derivatives are given, respectively, by 

K,, = ^^ = ^ [hab + daN'^ Kb + dbN' hac + Och^b) , (9) 
Kab = \ ^ = \ {kb + daN' hcb + dbN^ hac + dcKb) , (10) 
and are called the extrinsic curvatures. Note that 

Kab-Kab=^, (11) 

and 

Kab=^{VaVb + VbVa) , (12) 
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where Va = e^w* and is the covariant derivative with respect to hat [5]. 

As was stated above, the set of FPDs forms a gauge symmetry group of the system, 
and expresses a redundancy in the description of the system. This can be gauge-fixed 
arbitrarily according to convenience in describing the dynamics of a given system. Two of 
the useful gauge fixings are the foUowings: 

[ Comoving frame ] This is useful in describing elastic materials and is defined by the 
condition iV''(^,t) = 0, which leads to v'{^,t) = X'{^,t). 

[ Laboratory frame ] This is useful in describing fluids and is defined by the condition 
X°-{C,,t) = from which follow the equations = A^'^, hab = Sab, P = PoV^, and 
_ _|_ ybQ^^a^ rjij^^ j^^^ cquatiou shows that the covariant derivative in this frame 
becomes the material (Lagrangian) derivative. 

4 Fundamental equations 

We now write down a set of equations which determine the time evolution of twelve variables 
X* (three), hab (six) and A^" (three) up to FPDs (three). Our strategy is to first write down 
equations which hold in both of the elastic and fluid limits (this can be done in specific 
frames), and then to make the equations covariant under FPDs. 

4.1 Equation of motion for hab{^,t) 

We first consider the limiting case of elasticity in the comoving frame (A^" = 0). We then 
have hab = because the natural shape does not evolve in time for elastic materials. Its 
FPD-covariant expression is Kab = 0, and thus we expect that nonvanishing Kab represents 
genuinely plastic (i.e. nonelastic) deformations. However, one can show that its trace part 
K = h""^ Kab always vanishes due to the mass conservation [5] : 

i? = 0. (1) 

In fact, in the laboratory frame, K is expressed as (l/p)(p + 9a(pf°)) [5] and thus must 
vanish. 

The above consideration shows that only the traceless (shear) component of the extrin- 
sic curvature Kab expresses the changing rate of the bonding structure, which should be 
proportional to shear strains in the following form [5] : 
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Note that the second term on the left-hand side can actually be omitted because K = 
0. The inverse of the proportional coefficient, r, has the dimension of time, and should 
correspond to the relaxation time. 

In order to check the consistency of this equation (named the rheology equation in [5]), 
we take two extreme limits. In the elastic limit (r — oo), we obtain Rat = 0, which again 
means that there occurs no time evolution of intrinsic metric. In contrast, in the fluid limit 
(r —7- 0), we have Eab oc hat, or hat oc hab- This means that the material acquires only 
homogeneous compression around each point, which is actually the characteristic property 
defining fiuid. 

4.2 Equation of motion for 

The dynamics of should be expressed as Euler's equation 

paa = -V'na. (3) 

The leading form of the stress tensor Tab in the derivative expansion can be determined by 
the following requirements: 

• Tab be symmetric and covariant under FPDs. 

• Tab be linear in the strain Eab and the spatial derivative of the velocity, VaVb- 

These requirements imply that Tab is a linear combination of the irreducible components 
of Eab and Kab = {l/2){VaVb + VbVa) (scc (^^): 

Tab = -2/x - - £ hab - 27 K, - C Khab , (4) 

K 

where 

E = h'''Eab, e;^,=Eab-lehab, K = h'''Kab, = Kab - ^ Khab ■ (5) 

The set of equations ([I]), ([2]), ([3]) and (jlj) give nine independent equations, and thus 
can describe the time evolution of twelve dynamical variables up to FPDs. In order to 
see the meaning of the coefficients in (jl]), we consider two extreme limits of elasticity and 
fluidity. In the elastic limit, we have Kab = and thus K'^^^ = {^DEab/ Dt){l + 0(5)), 
K = {^De/ Dt) (1 + 0(e)). Then the stress tensor becomes 

Tab = -2/x EU--ehab-2^^-C ^hab + 0{e') , (6) 
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which shows that the four terms represent, respectively, shear and bulk moduli, and shear 
and bulk frictions. On the other hand, the fluid limit is realized by considering at time 
scales much longer than the relaxation time r. By using the rheology equation in the form 
= rKh{^ + 0(e)), Eq. © becomes 

Tab = -2v Kb - C Khab --ehab + 0{e^) , (7) 

K 

where 77 = 7 + /ir and ( are shear and bulk viscosities. The last term expresses the pressure 
with K being the bulk compressibility. 



5 Dispersion relations 

In order to get dispersion relations for the set of equations ([I])-(j4]), we expand solutions 
around the hydrostatic solution in the laboratory frame = hab = Sab) as 

hat = 6ab + Ahat, N'' = + AN\ (l) 

We consider plane waves propagating in the X^{= ^^)-direction: 

Ahabi^, t) = Ahabik, a;)e^^-«'-^-* , AiV'^(e, t) = AN''{k, a;)e^^«'-^'^* . (2) 

The equations can be separated according to the helicity, and we have the following equa- 
tions: 

• hehcity (shear modes) {AA^^ ± iAN"^ , Ah-^i =f iA/iga}: 



+ i 



1 1 4 



1 



Po 



1 / . 4 \ /I 4 



- + -{C + :,l)k' uj' U + -^] + [- + -p) k'co~^k' = (3) 



r V 3 / Vk 3 



PqKT 



T po\ 3 

fluid limit (r ^ 0) : + —(( + -7) k^u - — P = (4) 
Po ^ 3 / PoK 

% ( 4\ 1/1 4\ 

elastic limit (r -> 00) : u? ^ ( C + -7 ) ( - + -/i I A;^ = (5) 

Pn\ 3/ PnVfi; 3/ 



Pq\ 6 / Pq\K 

helicity 1 (sound modes) {AA^^ , A/i33 , Alin + A/i22} 



+ i 



^- + ^e)u-^^-±^e = o (6) 

r po ^ Po T- 



fluid limit (r -> 0) : ioj - —k^ = (7) 
Po 

elastic limit (r 00) : + i — k'^tu k'^ = (8) 



Po Po 



helicity 2 (scalar modes) {Ahu — A/i22 =F 2^ A/112}: 



^ = ~ (9) 
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6 Conclusion and discussion 



In this work we have shown that viscoelasticity can have a universal description by in- 
troducing an intrinsic metric to express bonding structures and by adopting FPDs as an 
underlying gauge symmetry. We have seen that FPDs naturally interpolate between elas- 
ticity and fluidity, and that covariance under FPDs uniquely determines the dynamics in 
the first order of strains. 

The next step should be to extend our framework to the systems with more than one 
relaxation time and with multiple components, in order to describe more realistic materials. 

We conclude with a comment that in the sound modes the diffusion equation becomes 
a telegraph- like equation when the relaxation time is taken finite (see Eq. (^). This 
means that information propagates only inside the "light cone", and thus the causality 
problem (inherent in diffusion equations) disappears when the effect of elasticity is taken 
into account. This property should be taken over also in the "relativistic extension of 
viscoelasticity" [B], and would serve as a causal completion of relativistic fluid mechanics. 

Acknowledgments 

MF thanks the organizers of the QTS6 conference, especially Sumit Das and Al Shapere, for 
inviting him to give a talk. Elucidating discussions with Michael Berry and Yuho Sakatani 
are also appreciated. This work was supported by the Grant-in- Aid for the Global COE 
program "The Next Generation of Physics, Spun from Universality and Emergence" from 
the Ministry of Education, Culture, Sports, Science and Technology (MEXT) of Japan. 
This work was also supported by the Japan Society for the Promotion of Science (JSPS) 
(No. 20-892) and by MEXT (No. 19540288 and No. 18540264). 

References 

[1] R. M. Christensen, Theory of Viscoelasticity, Academic Press (1971). 

[2] E. C. Bingham, Fluidity and Plasticity, Kessinger Publishing (2007). 

[3] L. D. Landau and E. M. Lifshitz, Fluid Mechanics, Butterworth-Heinemann (1987). 

[4] L. D. Landau and E. M. Lifshitz, Theory of Elasticity, Butterworth-Heinemann (1986). 

[5] T. Azeyanagi, M. Fukuma, H. Kawai and K. Yoshida, Phys. Lett. B 681 (2009) 290 
|arXiv:0907.0656l [hep-th]]. 

[6] M. Fukuma and Y. Sakatani, in progress. 



10 



